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While measuring the ratio Rk — r(^^j^|'^j)/r(A'*^^^^), the KLOE Collaboration has studied the 



radiative process Ke2y- The ratio of the width for the Ke2y decay with a positively polarized photon 
from structure-dependent radiation to the inclusive K^2{y) width is found to be 1.484(68) x 10^^. 
The observed radiation spectrum agrees with predictions from chiral perturbation theory and is in 
contrast with predictions based on the light front quark model. This result reduces the contribution 
to systematic uncertainties on measurements of Rk- In a separate study, KLOE has measured the 
ratio of the radiative Kej,y decay width to the inclusive K^^iy) width to be 924(28) x 10^^. The 
distribution in energy and angle of the radiative photon has been analyzed in an attempt to isolate 
the signature from interference of the inner-bremsstrahlung and structure-dependent amplitudes. 
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1. K^^e^v{v)y{Ke2y) 

The Kf.2 decay is strongly helicity suppressed. Its rate is therefore a sensitive probe for minute 
contributions from physics beyond the Standard Model (SM). This is particularly true of the ratio 
Rk = r(A'g2(y))/r(A'^2(7))» which can be calculated in the SM without uncertainties from strong- 
interaction dynamics. A recent calculation, which includes i^{e^p'^) corrections in chiral pertur- 
bation theory (ChPT), gives Rf^ = 1A11{\) x 10"^ |jl]]. Deviations of of up to a few percent 
are possible in minimal supersymmetric extensions of the SM with non- vanishing e-1 scalar-lepton 
mixing [|]. 

Additional photons in the Kg2 final state can be produced via internal bremsstrahlung (IB) or 
structure-dependent (direct) emission (SD); the spectrum of the latter is sensitive to the structure 
of the kaon. Interference between the IB and SD processes is negligible [^]. By definition, Rk 
is IB inclusive and SD exclusive. However, since events with IB photons cannot be distinguished 
from those with SD photons, to compare data with the SM prediction at the percent level or better, 
sufficiently precise knowledge of the SD component is required. 

Because of the helicity suppression of the IB channel, the contribution from the SD channel 
to the total width is approximately equal to that from IB Previous measurements give the 
SD rate with a relative uncertainty of 15% [Q]. The KLOE Collaboration has recently submitted 
for publication a measurement of Rk to within 1.3% ^ ^, obtaining Rk = 2.493(31) x 10^^. 
The study of the SD radiation in Ke2y decay presented here was conducted in order to reduce 
the contribution to the systematic error on the KLOE measurement of Rk from the uncertainty in 
the SD rate, from 0.5% to 0.2%. The NA62 Collaboration has recently announced a preliminary 
measurement of Rk with a 0.6% error [0]. The uncertainty in the SD rate is the second-largest 
contribution to the overall systematic uncertainty on this result, so that the present study may help 
to reduce the uncertainty on the NA62 value of Rk as well. 

Decays with SD radiation can proceed through vector and axial transitions, with effective 
couplings V and A, respectively: 

my " 647^2 [(^+^) /sD+(x,3') + (^-A) /sD-ix,}-)] . (1.1) 

Here, x = 2E*/Mk and y = 2E*/Mk are the normaUzed photon and electron energies in the kaon 
rest frame; both lie between and 1. The form factors /sd+ and /sd- describe the hadronic- 
structure-dependent effects in the contributions to the decay intensity from the channels with posi- 
tive and negative photon polarization, respectively. In the SD+ channel, the photon is preferentially 
emitted in the direction opposite to that of the e^, and the momentum spectrum is hard, peaking 
above the endpoint of the p* spectrum in K^^ decays at 228 MeV. {p^ is the momentum in the 
decay rest frame.) In the SD — channel, the photon is preferentially emitted in the direction parallel 
to that of the e^, and the has a softer momentum spectrum. This is illustrated in the Dalitz 
plots for each of the two contributions in Fig. |l|. Because the contribution from the SD— channel 
is completely submersed in ^e3(y) background, the present measurement is exclusively focused on 
the SD+ contribution, and events with p* < 200 MeV are eliminated to reduce background. 

The starting sample of Kg2[y) events (inclusive with respect to the presence of additional pho- 
tons in the final state) is the same as that used in the KLOE measurement of Rk [0, §]■ Single-prong 
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Figure 1: Plots of expected decay intensity in the {E*,p*) plane corresponding to the contributions from the 
SD+ (left) and SD- (right) channels, evaluated from the &'{p'^) ChPT calculation of Ref. [|]. 

vertices in the drift chamber are examined if the inbound track originates from the interaction point 
and has momentum compatible with the two-body decay momentum in ^ ^ K^K (127 MeV), 
and if the outbound track has pi^p > 1 80 MeV. The mass of the outbound track (the lepton candi- 
date) is then reconstructed as Mf^^ = {Ek — IPmissI)^ ~ /'kp' ^^^^ the kaon energy and Pmiss the 
missing momentum at the vertex. Stringent track quality cuts on both the inbound and outbound 
tracks increase the ratio of ^^2(7) events to background events from 1/1000 to 1/20. For the anal- 
ysis of Rx, the lepton-candidate track is extrapolated to the KLOE electromagnetic calorimeter, 
and a neural network is used to combine E/p, time-of-flight, and longitudinal energy deposition 
information for the purposes of n/e separation. 

For the analysis of Ke2Y, a hard cut is made on the neural-network output. The radiative photon 
must be explicitly detected as a cluster not associated to any track, with f™' > 20 MeV as measured 
in the calorimeter. This reduces background from activity in accidental coincidence with the event 
and from cluster fragments. Furthermore, the times of arrival of the lepton and photon clusters 
must be consistent with each other (given the event topology) to within 2a, or about 100 ps. This 
rejects K — > tttt*' events, in which the "lepton" is actually a pion, and travels with j8 = 0.8 instead 
of j8 1, so that its cluster arrives late. It also rejects events with accidental clusters, since the 
arrival times of such clusters are distributed evenly. 

After the above cuts, the dominant backgrounds are from K^2{r) events for M^^ < 20000 MeV^, 
and from K,^y) events for > 20000 MeV^. Ke2Y events with /?* < 200 MeV are submerged 
in background, as suggested by the left panel of Fig. ^ For the purposes of this analysis, the 
signal sample of SD+ events is fiducially defined as the sample of identified Ke2Y events with 
E* > 10 MeV, cos d* < 0.9 (d* is the angle of the photon momentum with respect to that of the 
e^), and p* > 200 MeV. Simulations indicate that this selection has a 90% acceptance for true SD+ 
events, a 2% acceptance for SD— events, and contains a 1% background of IB events. 

An additional variable that is useful in isolating signal events is AEy, the difference between 
the energy of the radiative photon as determined by kinematic closure using the tracks in the event, 
and that from the direct measurement of the energy with the calorimeter. On average, the photon 
energy is determined with a resolution of about 12 MeV in the former case and about 30 MeV in 
the latter. Since the Ke2y kinematics are built into this variable, it is particularly useful in separating 
Ke2Y from Kg^(^yj events, as seen in the right panel of Fig. |^ 
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Figure 2: Left: Distributions of Mj^p for K^iy candidates and backgrounds. The distributions for simulated 
Ke2y events with E* > 10 MeV and p* > 200 MeV (p* < 200 MeV) are shown in red (blue). The expected 
contribution from A'g3(y) events is shown in green; that from all sources combined is shown in black. The 
peak at M^,^ « 10000 MeV^ w is from events. The points with error bars are the data. Right: AE^/ a 
for simulated Ke2y events (yellow shaded histogram) and /T^sfyj events (open histogram). 
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Figure 3: ArsD+ (^^c27) /r(A'^2(y) ) for tiins in E*, with model predictions. Left: Overlay (not a fit) 
of &{p'^) ChPT prediction. Middle: Fit with ff{p^) ChPT prediction to obtain normalization and slope 
of vector form factor Right: Overlay (not a fit) of LFQM prediction with parameters as in Ref. Red 
histograms show expected contribution from IB. 



The number of Ke2y events in the SD+ sample is determined via log-Iikelihood fits to the 
distribution of candidate events in the (Mj^p, AEy) plane in five different bins in E*. The first bin 
in E* extends from 10 to 50 MeV; the remaining four adjacent bins are 50 MeV wide. The fit 
function is a linear combination of Monte Carlo (MC) distributions for Ke2Y events with kinematics 
as generated satisfying the SD+ selection (signal events), K^iy events as generated not satisfying 
the SD+ selection (dominantly IB events), K^2{y) events (including those with background photons 
in accidental coincidence), and events from other background sources (mainly K^-^i^y^). All events 
in the input distributions satisfy the analysis cuts for the SD+ selection as reconstructed. The free 
parameters of the fit are the weights for each distribution. 

Figure ^ summarizes the results of this analysis. The points in each panel show the measured 
values of ArsD+(^<s'27)> normalized to r(A'^2(7))» for each of the five bins in E*. Summed over 
all bins in£;,A^sD+(W = 1378 ±63, which gives rsD+(^.2r)/r(^^2(r)) = 1.484(66)st(16)sy x 
10^^. This agrees with the prediction from the KLOE MC, 1.447 x 10^^, which is based on 
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the &{p'') ChPT calculation of Ref. [|], in which V and A in Eq. (jTl]) are the constant values 
0.095/mjf and 0.043/mif . In the left panel of the figure, the ^{p^) result is overlaid with the data; 
as an index of the agreement, ;^^/ndf = 5.4/5. On the basis of this agreement, the KLOE MC is 
determined to be accurate in its description of the SD radiation in Kg2y decay to within the 4.6% 
relative uncertainty of the measurement. This reduces the corresponding contribution to the relative 
uncertainty on the value of Rk to 0.2% [||, ^. 



At 6{p \ the vector coupling V in Eq. ( pTTI ) depends linearly on the photon energy x, with 



V = Vo[l + X{l-x)], Vo = OMl/niK and A = 0.4; the axial coupling is constant, A ?a 0.034/mjf 
1^. The middle panel of Fig. || shows the data together with the results of a fit using the ff{p^) 
prediction. Since there is no sensitivity to the SD— component of the radiation, the quantity (V —A) 
in Eq. ( pTT] ) is fixed to its ^{p^) numerical value, and the fit parameters are (Vq +A) and A. The 
results (Vo +A) = 0.125(7)/m/f and A = 0.38 ±0.21 with p = -0.93 and ;t:^/ndf = 1.97/3 are 
obtained. The right panel of Fig. ^ shows the data together with the results of a fit incorporating 
the light front quark model (LFQM) which features a complicated dependence of V and A on 
X, with V=A = Oat;c = 0. The data are clearly in disagreement with the prediction of this model, 
as the overlay gives ;^^/ndf = 127/5. 

2. KL^n^e^v{v)y{Ke3y) 

The study of the Kg^j decay provides the opportunity to quantitatively test predictions from 
ChPT. In addition, precision measurements of the fully inclusive A'a(y) decay rates are needed 
for the determination of the CKM matrix element |V„s |. Such measurements require an accurate 
understanding of the radiation in these decays. 

As in the case of Kf.2j decays, the relevant kinematic variables for the study of Ke},y decays are 
E* and Qy. The IB amplitudes diverge for E* 0. The IB spectrum in 0* is peaked near zero as 
well, since m^. ?a 0. The IB and SD amplitudes interfere. The contribution to the width from IB-SD 
interference is 1% or less of the purely IB contribution; the purely SD contribution is negligible. 

In the &{p^) ChPT treatment of Ref. [^, SD radiation is characterized by eight amplitudes, 
{Wi, A,}, which in the one-loop approximation are real functions. These terms have similar photon 
energy spectra, with maxima around = 100 MeV. This suggests the following decomposition of 
the photon spectrum: 

The function f{E*) summarizes the effect of the SD amplitude; the parameter {X) measures its 
strength. The ^{p^) ChPT estimate is {X) = — 1.2ib0.4. The low-energy constants (LECs) for the 
^{p^) terms are unknown. An educated guess of their size leads to the assignment of an uncertainty 
on {X) of 30% of the ^(/) result. 

KLOE has studied the double differential rate in Kej,y decays, d^T/dE* ddy, and measured the 
ratio R, conventionally defined as 

T{Keiy\E* > 30 MeV, 6; > 20°) 
R= ^ ^ ^ ^. -, (2.2) 

r(^e3(7)) 
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The value of this ratio has been computed at ^{p^) in ChPT, leading to the prediction [10] 

/? = (0.963 + 0.006 (X)± 0.010) x 10"^ (2.3) 

For {X) = —1.2, R = (0.96 ±0.01) x 10^^, as quoted in Ref. [^. The simultaneous measurement 
of R and {X) allows a precise comparison with the theory, in large part avoiding complications 
from the uncertainties on the LECs for G{p^\ 

A first attempt to measure both R and the magnitude of the SD contribution was performed by 
the KTeV collaboration in 2001 [jn]]. They obtained R = (0.908 iO.OOS^o oij) x lO"^. the analysis 
of the SD radiation was complicated by the use of a cumbersome theoretical framework. The 
KTeV data were subsequently reanalyzed using more restrictive cuts that provide better control over 
systematic effects, but which reduce the statistics by a factor of three. The more recent KTeV result 
is /? = (0.916 ± 0.017) x 10^2 No further attempt was made to isolate the SD component. In 
2005, NA48 measured R = (0.964 ibO.OOS+oJJog) x 10"^ ^gain without attempting to isolate 
the SD radiation. 

The KLOE MC generates only radiation from IB, so the i&{p^) generator of Ref. [P] was used 



for the present analysis. The KLOE generator ||14|], makes use of a resummation in the soft-photon 
limit to all orders in a of the ^{p^) amplitude for single-photon emission. It is claimed to describe 
the IB photon spectrum at the level of ^\%. While this is appropriate for inclusive decay -rate 
measurements at the 0.1% level, the SD contribution is about 1% of the IB contribution, so that the 
accuracy level of the KLOE IB generator is of about the same order as the SD contribution itself. 
Therefore, in this analysis, the generator of Ref. [^] is used to obtain the photon spectrum from IB 
as well as from SD. 

The analysis is fully described in Ref. [|l5|]. The criteria used to select an inclusive sample of 
Ke?,Y events are the same described in Ref. [jl^]. Vertices in the drift chamber along an expected 
hne of flight reconstructed from a — > tt+tt^ decay are identified as candidate Ki decays. Loose 
kinematic cuts are applied to remove background from Ki n^n^n^ and Kl ^ n^n^. A large 
amount of ^^3(y) background is rejected using the variable Aj^^, the lesser value of l^miss — /'miss I 
calculated in the njJ. and /itt track mass hypotheses. Events are retained only if Aj^^ > 10 MeV 
The difference between the times of flight of the hypothetical and tt^ tracks is used to purify 
the sample to about 99.3%. 

Signal KeSy events are selected from within the inclusive sample. Events are first required to 
contain a calorimeter cluster not associated with any track. The arrival time of this cluster must be 
within 8 a of the time expected on the basis of the position of the decay vertex and the photon time 
of flight. If there is more than one photon candidate, that with the most compatible timing is used. 
Once the photon has been identified, the requirement 0* > 20° is enforced. The photon energy in 
the laboratory system is calculated from the track momenta and the photon cluster position, using 
the constraint that /My = 0. The resolution obtained is '^1 MeV, or about a factor of ten better than 
that of the energy measurement from the calorimeter. To check the agreement between MC and 
data on the reconstruction performance for quantities used for photon selection and to correct the 
MC efficiency, Kl — > n^n^n^ events are used as a control sample. 

Neural-network discriminators are used to reduce backgrounds from Ki n^n^n^ and ^^3(y). 
The neural network used to identify Kl n^n^n^ events makes use of E*, Q*, the track momenta. 
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Figure 4: Results of fit to {Ey, By) distribution: dots show data, dark gray region shows contribution from 
Ke2,y events (IB and SD) satisfying kinematic cuts, white region shows contribution from K^^y events not 
satisfying cuts, light gray region shows contribution from background. 



the missing momentum, and M^, the invariant mass of the photon-neutrino pair. The neural net- 
work used to identify ^^3(y) events is based on the track momenta, the calorimeter energy mea- 
surement, and the depth of penetration into the calorimeter of the cluster centroid. 

To count the numbers of K^^y events from IB and SD, a fit to the experimental distribution 
in {Ey,dy) is performed using the sum of four independently normalized MC distributions: the 
distribution for Ke3y events from IB that satisfy the cuts on E* and 6* as generated; the distribution 
corresponding to the function /{£*) in the second term of Eq. ( |2.1| ), representing the modification 
of the spectrum from SD events; KeSy events from IB that satisfy the kinematic cuts only as recon- 
structed (and not as generated); and physical background from Kl — > n^n^'nP and ^^3(7) events. 
The free parameters of the fit are the number of IB events, the effective number of SD events (the 
integral of the spectral distortion induced by the IB-SD interference), and the number of Ke-iy events 
not satisfying the kinematic cuts. The magnitude of the background contribution is fixed using the 
MC. Figure ^ shows the result of the fit. The two-dimensional distributions are plotted on a single 
axis; the E* distributions for each of the eight slices in Q* are arrayed sequentially. The fit gives 
;t: Vndf = 60/69 {P = 11%), with Mb (^.37) = 9083 ±213, A^sD(^rfr) = - 102 ± 59, and an addi- 
tional 6726 ± 194 Ke'iy feed-down events not satisfying the kinematic cuts. The negative value for 
the effective number of SD events is a result of the destructive interference between the IB and SD 
amplitudes. The value R = 924(23)st(16)sy x 10^^ is obtained. The value of the parameter {X) 



defined in Eq. (gjj) is derived from the result of the fit, taking into account the small difference in 
the overall detection efficiencies for IB and SD events. The result is (X) = —2.3 ± 1.3st ± l-4sy. 
The la contour is illustrated as the open ellipse in Fig. S. 



The dependence of R on {X) of Eq. ( |2.3| ) is shown in Fig. g as the diagonal shaded band. This 
dependence can be used to further constrain the possible values of R and (X) from this measure- 
ment, giving the 1 a contour illustrated as the filled ellipse in the figure. The constraint is applied 
via a fit, which gives R = 944(14) x 10"^ and (X) = -2.8 ± 1.8, with correlation p = 72% and 
;^^/ndf = 0.64/1 {P = 42%). This result represents an improved test of ChPT with respect to that 



obtained using the measurements of R from Refs. [12, |13| 



Finally, to test the accuracy of the KLOE IB generator [14], fits to the data with no SD com 



7 



Recent KLOE results on radiative kaon decays 



presented by Matthew Moulson 



0.98 



0.96 



0.94 



0.92 



0.90 









KLOE 

ChPTconsujta / 




, <- KLOE 











-10 



Figure 5: KLOE la contours in the (/?, (X)) plane from fit to the (Ey, By) distribution (open ellipse), and 
same results when combined with constraint from ChPT (filled ellipse). Results from KTeV ^ and NA48 
[ pj] ] are also shown, as well as the dependence of Ron {X) used as the constraint. 



ponent have been performed. The results R = 921(23)st x 10~^ for the KLOE generator [14] and 
R = 925(23)st X 10"^ for the ^{p^) generator of Ref. [Q] have been obtained. This confirms the 
reliability of the KLOE generator for IB events. 
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